Self accelerating electron Airy beams 
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We report the first experimental generation and observation of Airy beams of free electrons. The elec- 
tron Airy beams are generated by diffraction of electrons through a nanoscale hologram, that imprints a 
cubic phase modulation on the beams' transverse plane. We observed the spatial evolution dynamics of an 
arc-shaped, self accelerating and shape preserving electron Airy beams. We directly observed the ability of 
electrons to self-heal, restoring their original shape after passing an obstacle. This electromagnetic method 
opens up new avenues for steering electrons, like their photonic counterparts, since their wave packets can 
be imprinted with arbitrary shapes or trajectories. Furthermore, these beams can be easily manipulated us- 
ing magnetic or electric potentials. It is also possible to efficiently self mix narrow beams having opposite 
signs of acceleration, hence obtaining a new type of electron interferometer. 
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An arc of light is a caustic phenomenon 1 1] abundant in 
nature; ranging from the rainbow in the sky|2], to the bright 
light patterns that appear at the sea floor when sun shines 
on the rippling waves of water|l]. In the framework of 
quantum mechanics, Berry and Balazs|3] found a unique 
wave packet of a massive particle in the form of the Airy 
function|2]; its counter-intuitive properties are revealed as 
it propagates in time or space. It preserves its shape de- 
spite dispersion or diffraction and continuously self accel- 
erates, with no force applied, along an arc-shaped trajec- 
tory. Nearly 30 years later, this wave packet known as 'Airy 
beam' was realized by Christodoulides|4] and co-workers 
in the optical domain; later it was generalized to accelerat- 
ing optical beams with an arbitrary spatial shape| 5] and any 
curved convex trajectory! 6]. Here we report the first experi- 
mental generation and observation of the Airy beams of free 
electrons. The electron Airy beam is generated by diffrac- 
tion of electrons through a nanoscale hologram|7, 8J, that 
imprints a cubic phase modulation on the beam's transverse 
plane. We observed the spatial evolution dynamics of an 
arc-shaped, self accelerating and shape preserving electron 
Airy beams. We directly observed the ability of electrons 
to self-heal restoring their original shape after passing 
an obstacle. This electromagnetic method opens up new 
avenues for steering electrons, like their photonic coun- 
terparts, since their wave packets can be imprinted with 
arbitrary shapes^ or trajectories 1 6]. Furthermore, these 
beams can be easily manipulated using magnetic or electric 
potentials [n], [n|. It is also possible to efficiently self mix 
narrow beams having opposite signs of acceleration, hence 
obtaining a new type of electron interferometer. 

When a particle's wave packet (representing its probabil- 
ity density), is given the initial shape of an Airy function, 
it does not spread and continuously self accelerate as it 
evolves in time. This was theoretically predicted by Berry 
and Balazs in 1979 [3] as they showed that the Airy wave 
packet is a solution to the Schrodinger equation for a free 
particle. The evolution of a slowly varying amplitude light 




FIG. 1. Electromagnetic generation of an electron Airy beam, (a) 

An electron beam is transmitted through a nanofabricated holo- 
gram, with a cubic phase modulation. It is then focused by a mag- 
netic lens. An electron Airy wave packet is formed at the back focal 
plane and recorded as it evolves. An electron Airy beam is a narrow, 
shape preserving, self accelerating beam. It also self-heals after 
encountering obstacles, (b, c) TEM micrographs of the nanoscale 
holograms (b) 2D Airy on a spatial carrier frequency, (c) 2D Airy 
without a carrier, (d, e) Experimental wave packet micrographs of 
2D and ID electron Airy beams. 



beam in space is analogous to the evolution of the wave 
packet of a massive particle in time; much as the paraxial 
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Helmholtz equation resembles the SchrQdinger equation. 
This analogy led Christodoulides and co-workers 1 4] to dis- 
cover and experimentally realize the 'accelerating Airy op- 
tical beams', the optical analog of the quantum Airy wave 
packet of a massive particle|3j]. The accelerating Airy wave 
packet does not contradict the Ehrenfest theorem as its 
center of mass propagates in straight lines, but the high- 
est probability amplitude of the particle's location, i.e., the 
strongest lobe of the Airy function preserves its shape and 
stays localized around a parabolic trajectory in space, sim- 
ilar to that of a freely propagating projectile experiencing 
a transverse accelerating force. This intriguing propaga- 
tion dynamics is a caustic wave phenomena which can be 
understood by ray analysis: multiple rays emerge sideways 
from the edge lobes of the initial wave packet area and coa- 
lesce along a curved boundary 1 12]. Various applications fol- 
lowed the discovery of optical Airy beams including micro- 
particle manipulationlljj], generation of plasma channels 
in air and water [3, surface Airy plasmons fljUl^l and 
applications in lasers HE 13 and in nonlinear optics Eoll . 
However, all these applications relied on the wave-packet 
of zero mass photons. We report here the first generation 
and observation of Airy wave packet|3] of free electrons, or 
'electron Airy beams' (generated without using photons). 
We demonstrate a new electromagnetic technique to gener- 
ate electron Airy beams analogous to the optical method|4], 
as we exploit the recent advances in nanoscale hologram 



It is then clear from Eq. [2] that |*| 2 preserves its shape and 
accelerates in the transverse coordinates. The acceleration 
is described by|H x(z) = z 2 /4fc|x 3 , y(z) = z 2 /4A;|y 3 . The 
aforementioned ideal Airy beam carries an infinite amount 
of energy, wheres the feasible Airy beam is truncated, hav- 
ing a finite energy. The finite Airy beam is obtained by mul- 
tiplying the Airy function with an exponential or Gaussian 
window|4]. Nevertheless, over a finite distance, the finite 
Airy beam exhibits all the special characteristics of the infi- 
nite Airy beam such as low diffraction, free acceleration and 
self-healing. 

In optics, finite Airy beams were experimentally obtained 
by passing a Gaussian beam through a phase mask imprint- 
ing a cubic phase modulation in the transverse direction|4], 
followed by an optical Fourier transform (notice that the 
Fourier transform of a function having a cubic phase modu- 
lation results in an Airy function) . In our experiment we uti- 
lize a field emission gun transmission electron microscope 
(FEG-TEM) , operating at 200 kV. The de-Broglie wavelength 
(including relativistic correction) in this case is approxi- 
mately 2.5 pm. Therefore to view the inner structure of Airy 



fabrication techniques 0, Si- We supply a robust method to 
manipulate the electrons by trajectory, as their high peak 
intensities (indicating the probability density for locating 
the electron) self accelerate in a parabolic path. Further- 
more, this concept can be generalized because recent stud- 
ies have shown that optical Airy beams can be designed 
with arbitrary spatial shapes 1 5] or propagate along arbitrary 
trajectories [6]. 

The spatial evolution of the envelope * of an electron's 
wave function, can be expressed by the paraxial Helmholtz 
equation (see Supplementary Information Sec. 1): 

V^+2ifc B Jy * = , (1) 

where V\ — d jx 2 + d jy 2 is the transverse derivative and 
kB — p/ti = 2n/?iB is the de-Broglie wave number of the 
electron. This equation has the same form as that of the 
Schrodinger equation. However, rather than measuring the 
evolution of the envelope of the electron wave function in 
time, we measure it in space. This is another manifesta- 
tion of the analogy, widely used in optics, between beam 
diffraction in space and a pulse dispersion in time. When 
the initial probabilistic wave function of the electron is an 
Airy function (Ai), *(x,y,z = 0) = Ai(x/x ) • Ai(y/y ), where 
Xo,yo are characteristic length scales; the general solution 
to the wave packet is given 
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beams as well as the spatial separation of different orders, 
we used nanoscale holograms (see Supplementary informa- 
tion Sec. 2). The emitted electrons pass through the holo- 
grams which adds a transverse cubic phase: 
exp(f0(x,y)) = exp(/a x x 3 )-exp(/a y y 3 ), to the wave func- 
tion. Our hologram design method was to construct a bi- 
nary diffraction grating with the following shape I22LI23II: 

h{x,y) = ~ fo (sgn{cos[27rx/A + a x x 3 + a y y 3 ] + D cyc i e } + 1). 

(3) 

In this way, a cubic phase is imposed on a carrier frequency. 
The carrier period is A, ho is the ridge height of the binary 
phase mask and D cyc i e is an arbitrary duty cycle factor. The 
micrographs of these nanoscale holograms are shown in 
Fig. QJ>, c. The most important holograms for our exper- 
iment are: A 2D Airy lattice structure with carrier (AWC): 
A x = 400nm, a x , y = 2- 10 10 -27r/(400- lO" 9 ) [1/m 3 ] shown 
in Fig. QJ). A simple periodic Bragg lattice (BR) structure 
with a period A^ = 100 nm. Fourier transform of the mod- 
ulated wave function is done by using a set of magnetic 
lenses, so that an electron Airy beam is obtained at the back 
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FIG. 2. Acceleration of Airy beams, (a) The measured trajectories 
of Airy orders (1,-1) (blue curves) fits the simulation (gray back- 
ground). Airy orders accelerate anomalously inward (to the center 
of the axis). However the carrier (black curve) deflects in the usual 
manner outward with a constant transverse velocity, (b-e) Micro- 
graphs revealing the evolution and acceleration of electron Airy 
beams v.s numerical simulation results (d) Airy orders can collide 
(in x axis) , (e) or pass each other. Notice that the peak intensities of 
the different orders rotate clockwise (e.g. for the -1 order), or anti 
clockwise (for the +1 order). 



focal plane of the FEG-TEM. This method is analogues to 
the method used to generate optical Airy beams [4], and the 
only difference is that we manipulate electrons, rather than 
photons. The experimental Airy profiles of 2D and ID Airy 
beams are shown in Fig. [TJi, e. Since the measurement plane 
is located at a fixed position in our FEG-TEM, we varied the 
focal lengths of the magnetic lenses in order to observe the 
formation and evolution of the electron Airy wave packet in 
space as illustrated in Fig. [1^. The experiment details are 
given in Supplementary information Sec. 3. 

In the first experiment we measured the accelera- 
tion of electron Airy beams (results displayed in Fig. 
|2). We recorded profiles using a relatively large area 
(100 jum X 100 /im). The diffraction pattern from an AWC 
includes several diffraction orders. The zero diffraction or- 
der is not an Airy beam, therefore it rapidly expands and 
becomes invisible at the measurement plane. However, all 
the higher order modes are accelerating and shape preserv- 
ing (hence spatially localized) Airy beams. The acceleration 
of these beams can be shown in the x-z plane or the trans- 
verse x-y plane as presented in Fig. |2k, (b-e) respectively. 
We simulated the evolution of the electrons, by using beam 
propagation numerical simulation 1 20]. The numerical re- 
sults are in good agreement with the experimental results. 



Detailed explanation of acceleration measurement is given 
in Supplementary Information Sec. 4. In particular, it can 
be seen that the Airy beams of order +1 and -1 accelerate 
toward each other, eventually colliding, as shown in Fig. |2fcl, 
e. This represents a new way of interfering electron beams: 
The electron beam is separated into two diffraction orders, 
but these beams re-merge owing to their opposite accelera- 
tion directions. The measured deflection is 1.7 mm over an 
effective distance of 70 m. Notice that the acceleration can 
also be viewed in a single plane, as the observed Airy or- 
ders reach different transverse heights (see Supplementary 
Information Sec. 5). 

In the second experiment, we compared between elec- 
trons diffracting from an Airy lattice to the electrons 
diffracting from a reference periodic Bragg grating. We sim- 
ulated the evolution of electrons (Fig. |3^) diffracted from 
a periodic Bragg lattice and an Airy lattice with the same 
carrier period. (Notice that in this case only, we used dif- 
ferent input beam area for Airy and Bragg lattices, in order 
to visualize the trajectories). As seen, the diffraction from 
a Bragg lattice is normally outward, however, the diffraction 
from an Airy lattice, is anomalously accelerating inward. We 
recorded profiles (Fig. |3p-d) using a relatively small area 
(10 [imXIO fim), letting the zero order to slowly spread out 
thereby enabling visualization of the ratio between the evo- 
lution of spreading Bragg peaks compared to the shape- 
preserving Airy lattice peaks. The zero order of Airy lattice is 
the only diffraction order, which is not imprinted with a cu- 
bic phase, thus, it spreads in the normal manner. The zero 
order appears as a large circular spot at the middle of the 
frame. We can also see diffraction orders 1, -1, 2, -2. Notice 
that the intensity of Bragg peaks with opposite signs is not 
symmetric as expected. This is probably because the sam- 
ple was slightly tilted and not perfectly perpendicular to the 
electron beam. As can be seen, at an effective distance of 
52 m and 100 m the diffraction patterns from BR (as well 
as from the zero order of AWC) spread and become huge, 
while the non-zero Airy peaks stay confined along a very 
large effective distance (even after an effective distance of 
100 m) to a very small scale area. This difference is empha- 
sized when calculating the full width half maxima (FWHM) 
of the diffracted orders. FWHM of Bragg orders (1,-1) are: 
(b), (c), (d)= 1125, 5785, 8125 (im, however, FWHM of Airy 
orders (1,-1) are: (b), (c), (d)= 102, 104, 110 /xm. We can 
compare the FWHM of the diffracted orders: The ratio be- 
tween the FWHM, i.e., FWHM(Bragg) /FWHM(Airy) is about: 
11, 58, 81 respectively Notice that the profiles resemble the 
profiles given in Fig. [2] but the Airy function in this case, is 
more truncated. 

We also measured the self healing properties of electron 
Airy beams, as presented in Fig. [4] For this purpose, we used 
a wire placed in the diffraction plane. The wire was conven- 
tionally used as a bi-prism in electron holography|24J], but 
in our case it was simply used to block parts of the beam. 
This experiment was conducted in the following way: in- 
creasing the current of the objective lens elevated the Airy 
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FIG. 3. Comparison between electrons diffracting from an Airy 
lattice to the electrons diffracting from a reference periodic 
Bragg grating, (a) The diffraction from a Bragg lattice is normally 
outward at angles|8]: a,„ = mA B /A. However, the diffraction from 
an Airy lattice, is anomalous because the lattice peaks accelerate 
inward (b-d) Experimental profile micrographs of different prop- 
agation planes, (notice that in both the Airy lattice and the Bragg 
grating the zero order looks the same). The Airy orders are very 
localized and maintain very high intensities compared with the 
Bragg orders, (e) An overlay of different profile micrographs. 



beam above the wire. Then we adjusted the wire to simul- 
taneously block the two Airy orders of -1 and 1. We blocked 
the different orders in a slightly different manner. We then 
gradually increased the current of the diffraction lens and 
observed their self-healing progress|9]. The wave packets 




FIG. 4. Self healing properties of electron Airy wave packets. Ex- 
perimental profile photograph versus numerical simulation. The 
diffracted orders 1 and -1 are blocked differently so their self- 
healing progress is different. The two orders are blocked with a 
glass wire. 



reconstructed their shapes after passing the blocking wire. 
As shown, the healing progress of the orders was different. 
Order 1 self-healed faster than order -1. We also simulated 
the self healing process and the numerical results are in a 
good agreement with the experimental results. 

We demonstrated experimentally, for the first time, non 
spreading, freely accelerating electron wave packets in the 
shape of an Airy function. This gives rise to a novel field 
of manipulating particles' trajectories by engineering their 
probability density wave functions. We experimentally 
demonstrated trajectories of electrons, accelerating in free 
space without induction of force. We observed shape pre- 
serving, highly localized electron Airy beams which didn't 
spread even after an effective length of 100 m. Such non- 
spreading Airy electron wave packets may possibly be use- 
ful for improving the resolution properties of TEM imag- 
ing, since they have an extremely large depth of focus and 
a resolution depth feature. Furthermore, we showed, for the 
first time, that these self accelerating electron beams can 
self heal and reform their original shape after passing an 
obstacle. We demonstrated an interesting feature of these 
beams, in which different diffraction orders with opposite 
signs of acceleration can be merged and possibly be used as 
a new type of electron wave interferometer. It may also be 
interesting to study with these beams electron spin interac- 
tion effects in the relativistic regime of the electron wave- 
packets, similarly to recent studies of electron vortices fllll. 
New possibilities may open up for manipulating and shap- 
ing the trajectories and the self healing properties of Airy 
beams, as electrons can be influenced by magnetic or elec- 
tric potentials | lOj] . 
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SUPPLEMENTARY INFORMATION 
1. Quasi-relativistic Schrodinger equation 

Ignoring spin effects, we may use the Klein-Gordon equa- 
tion instead of Dirac equation in order to find the wave 
function of an relativistic electron. For free space propaga- 
tion (without any external potential): 

32 

-h 2 —<S> = {mc 2 ) 2 <S>-c 2 % 2 V 2 <S>. (4) 
at 2 

Its electron wavefunction plane wave solution is: 

$(r,t) = $ ex'p{ip-r/h-iEt/ti) , (5) 

where E, p satisfy the dispersion relation 
E—^/c 2 p 2 —{mc 2 ) 2 . Classically, E — ymc 2 , P — fPmc, 
j = l/sjl-p* and p = v/c. 

For small angle diffraction p± <C p (paraxial approximation) 
we can look for a wave solution of the form: 

$(r,t) = <li{r L ,z)-exp{iPz/h-iEt/H) , (6) 

Assuming slowly varying envelope and neglecting the sec- 
ond derivative of z, from Eq. [4] we obtain the paraxial 
Helmholz equation: 

Vi + 2«fc B ^j* = , (7) 

where V\ — d jx 2 + d jy 2 is the transverse derivative and 
ka — plfl — jmv/Ti — 2ti/Xb is the de-Broglie wave number 
of the electron. This equation has the same form as the 
paraxial Helmholtz equation, which describes the propaga- 
tion of light beams in space (used to find the solution for 
Airy optical beams 1 4]), except that it includes the de-Broigle 
wave number. 

Setting z — fict results in: 

f h 2 2 d \ 

V 2 ,+in— * = 0, (8) 

y2jm e 1 3t J 

which is the same as the Schroedinger equation that was 
solved (in ID) by Berry and Balzas|3], with a relativistic cor- 
rection to the electron rest mass. 
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2. Nanoscale hologram preparation 

Our hologram design method was to construct a binary 
diffraction grating with the following shape Ei . E^l : 

h{x,y) = -h {sgn{cos[2nx/A + a x x 3 + a y y 3 ]+ Dcyde} + 1). 

(9) 

In this way, a cubic phase is imposed on a carrier frequency. 
The carrier period is A, h is the ridge height of the binary 
phase mask. D cycle is an arbitrary duty cycle factor. Set- 
ting the duty cycle factor < D cy cie < 1 to high values (0.8 
, 0.9...), gives a hologram with very narrow stripes. This 
method was chosen because it matches well with the high- 
est resolution writing mode of the focused ion beam (FIB). 
The nanoscale holograms were prepared by sputter deposit- 
ing 10 nm Au layer onto a 50 nm silicon membrane chip, 
then milling the desired pattern with a Raith IonLine fo- 
cused ion beam machine. The machine mills with a 35 keV 
Ga ions beam, and we used 2.6 pA ion beam current with 
step size of 2nm and dwell time of 1.54 msec. This gives an 
effective dose of 20000 pC/cm for single-pixel-line milling. 
These conditions milled through the entire gold layer and 
in addition about 20 nm of the silicon nitride membrane. 
The milling was done over areas of 30 fim X 30 fim as shown 
(only in the central part) in Fig. lb, c. This method enabled 
us to achieve high resolution writing, without breaking the 
membrane chip. This kind of resolution was necessary for 
our experiment. 

1. A 2D Airy lattice structure with carrier (AWC): 
A = 400 nm and a x , y = 2 • 10 10 • 2?r/(400 • lO" 9 ) [1/m 3 ] 
shown in Fig. lb. 

2. A 2D structure with no carrier (ANC): a x ,y = 
2 • 10 10 • 2tt/(400 • 10~ 9 ) with the same acceleration as 
AWC [1/m 3 ] as shown in Fig. lc. 

3. A ID structure with carrier: Ax = 400nm, and accel- 
eration only in y: a y = 2 • 10 10 • 27r/(400 • lO" 9 ) [1/m 3 ] 

4. A simple periodic Bragg lattice (BR) structure with a 
period Ax = 100 nm. It served as a reference scale for 
the measurement. 



3. Experiment details 

We used a Tecnai F20@ FEG-TEM and studied the gen- 
eration and evolution of electron Airy beams by varying 
the focal lengths of the TEM magnetic lenses 1 8]. We im- 
aged the electron wave packet as it evolved in different 
planes between the back focal plane and the image plane. 
The profiles were taken in FEG-TEM low angle diffraction 
mode (LAD), for which the objective lens is with low cur- 
rent (around 7 % of the maximal current). This enabled high 




FIG. 5. (color online). Micrograph of AWC, ANC BR. The afore- 
mentioned method of writing enabled us to achieve high resolu- 
tion, without breaking the membrane chip. This kind of resolution 
was necessary for our experiment. 

magnification of the back focal plane. Setting a high current 
in the magnetic condenser lens and using an area adjuster, 
gave a relatively low convergence electron beam. Then, by 
a 'Free Lens Control' software increasing the current in the 
diffraction magnetic lens only (which was the only lens we 
changed during the measurements), we imaged different 
planes following the back focal plane as shown in Fig. la. 
A Gatan camera model 694 captured the wave-packet evo- 
lution and the spatial propagation dynamics of the electron 
Airy beam. 

4. Acceleration measurement 

One of the challenges in this setup is to determine the ef- 
fective distance that the beam travels from the focal plane. 
We change this distance by decreasing the focal length of 
the magnetic lens each time. In order to calibrate the prop- 
agation (z) axis we made the following analysis: We mea- 
sured the Bragg diffraction patterns from the periodic refer- 
ence grating for different setting of the magnetic lens. The 
location of the first diffraction order is given by|8]: 



Peak x{BR) {z) = x s[BR) + - z. (10) 

This diffraction is measured using the camera, and since 
each pixel size is 25 microns, we can obtain the distance z 
in meters. As for the diffraction of the Airy wave through the 
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FIG. 6. (color online). For measuring the acceleration we used 
three structures : AWC, periodic Bragg reference (BR) and addi- 
tional structure only having an Airy modulation with no carrier 
[Airy with no carrier (ANC)]. We have also recorded the electron 
evolution from ANC and BR for same z (at the identical experimen- 
tal parameters) and analyzed the results.structure 



mask with carrier (AWC), it has both acceleration and linear 
velocity components in x direction, i.e. the peak location of 
the strongest lobe is given by lilll : 



Peak x ( Ai )(z)=x S ( A i) + 



A, 



-z + - 



(11) 



whereas in the Y direction we have only acceleration, hence 
the peak location is: 



Peaky(Ai){z) = 



1 



Therefore, the difference between then is: 



Peak x ( Ai) {z) - Peak y (Ai){z) = x s{M ) + 



Vx(Ai) 



(12) 



(13) 



which is also linearly increasing with a constant transverse 
velocity in z. Since the period here is A x (Ai) — 400 nm, vs. 
only A x (£j?) = 100 nm for the Bragg reference mask, the slope 
and the starting point x s in the case of the Airy beam should 
be 4 times lower. We can therefore verify the calibration by 
comparing the two measurements, as shown in Fig. [6}3. 

We can now compare the measured peak location in the X 
axis, when diffracted from the AWC mask with the theoret- 
ical prediction. The +1 and -1 diffraction orders are clearly 
seen, and they are accelerating toward the center. There is 
an excellent agreement with the theoretical prediction. In 
addition, we have the results from this mask in the Y direc- 
tion of orders +1 and -1. Since there isn't a carrier modu- 
lation in the mask's Y direction, the beam starts in this case 
from the axis and accelerates away from it, as shown in the 
red curve of Fig. [6^. We have also measured the diffraction 
patterns from the ANC mask. Here there is no carrier mod- 
ulation, neither in x nor in y, hence the beam starts from 



the axis and accelerates away from it (green line in Fig. [6} a. 
As expected, it nearly coincides with the results of the Y axis 
from the AWC mask. The measured deflection is 1 .7 mm ( 70 
pixels) over an effective distance of 70 m. 

5. The Airy lattice 



(a) Fourier plane Z=0 
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FIG. 7. (color online). Localized accelerating lattice. The accel- 
eration can be also viewed in one profile picture taken after the 
back focal plane. Different orders have different acceleration co- 
efficients, (a) Simulation of Fourier plane of an Airy lattice, (b) 
Simulation of a following plane. As can be seen each order ele- 
vate to a different height (in y). (c) Experimental profile showing 
this effect, (d) Airy orders have also different spatial distribution 
as shown with the different orders of the ID structure 

Currently, the acceleration measurement is done by tak- 
ing several profile pictures in different plains following the 
Fourier plane. The purpose of this section is to explain how 
the acceleration in the Airy lattice can be indicated by single 
plane (after Fourier plane) . 

The Airy lattice is a periodic lattice imposed with a cubic 
phase. A general formula to impose a lattice with an arbi- 
trary phase is (in binary structures) |3 : 



h(x,y) — -ho ■ sgn{cos[2nx/A + 0(x,y)]). 



(14) 



Any arbitrary phase can be imposed on this lattice. As 
an example the fork shaped structure 0, is generated 
by: <j)(x,y) — / c tan _1 (y/x), where l c is the topological 
charge. The Airy lattice is imposed with following phase: 
(p(x,y) — a x x 3 + a y y 3 . When electrons (or light) diffracts 
from the aforementioned binary structure, it decomposes 
to different orders; the complex amplitude of the m-th or- 
der diffracted beam is oc exp(/ • m(j)(x,y)). In the case of 
fork shape structure, series of vortices is obtained. Each 
vortex has a different orbital angular momentum, which 
equals to m ■ l c . The radius of each vortex increases as a 
function of orbital angular momentum, thus, the diffraction 
from a fork shaped structure generates a series of vortices 
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with different radii 0,0]. In the case of Airy lattice, each 
diffracted order is imposed with a different cubic phase: 
oc exp(; • m{a x x 3 + a y y 3 )), thus, each order experience dif- 
ferent acceleration. This can be viewed in a single plane as 
shown in Fig. [7]it was also observed experimentally Fig. [7J: 



This is a novel type of lattice. Although it is decomposed 
to different orders which propagate in different directions, 
each single order stays localized (unlike Bragg lattice) and 
anomalously accelerate. 



